Limitation of energy deposition in Classical N Body Dynamics. 
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Energy transfers in collisions between classical clusters are studied with Classical N Body Dynam- 
ics calculations for different entrance channels. It is shown that the energy per particle transferred 
to thermalised classical clusters does not exceed the energy of the least bound particle in the cluster 
in its "ground state" . This limitation is observed during the whole time of the collision, except for 
the heaviest system. 
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The question of energy deposition in nuclei during 
nucleus-nucleus collisions is of great importance for nu- 
clear matter studies. The maximum amount of energy 
which can be stored in hot equilibrated nuclei has been 
studied both experimentally [H Q C3- Gl 0- El and theo- 
retically 0, @ • Such studies have mainly been motivated 
by the determination of a plateau in the so called caloric 
curve (the evolution of the temperature with the excita- 
tion energy) which could be a signature of a first order 
liquid-gas phase transition 0, . 

Experimentally, a lot of work has been done to de- 
termine the amount of thermal energy stored in nuclei. 
For central collisions, large energy dep osits up to 22.5 
A.MeV have been determined 0, uM- Other studies 
have shown that the excitation energy in primary prod- 
ucts in multifragmenting systems is around 3-4 A.MeV 
QL 13 > for below the total available energy in the 
center of mass frame. This energy does not seem to evolve 
strongly with the incident energy. These two measure- 
ments seem to be in contradiction. Possible limitations 
of energy deposition could result from prompt emission 
of energetic ligh t charged particles at early times in the 
reaction [II El EH. 

Theoretically, the maximum energy that a equilibrated 
nucleus can withstand corresponds to the energy (or tem- 
perature) at which the surface tension vanishes. This is 
often characterised by a critical temperature T c whose 
value is around 16 MeV [7], |8|]. This temperature is 
linked to the equation of state of nuclear matter. The 
main drawback of such studies is that they assume that 
the system is fully equilibrated and hence do not take 
into account possible limitations coming from the reac- 
tion mechanism. 

The aim of the present article is to study energy deposi- 
tion during collisions between finite size systems in a well 
controlled framework. Results from the Classical N Body 
Dynamics code E3 wn ^ be shown and the mechanism of 
energy deposition in classical clusters will be studied. We 
will be interested in this paper by the energy transfered 
to thermally equilibrated clusters which corresponds to 
the energy that long-lived clusters can withstand. The 
paper is organised as follows. In a first section, the Clas- 
sical N Body Dynamics will be briefly described. The 
excitation energy in clusters will be shown for various 



systems and incident energies in the second section. The 
third section will be devoted to the mechanism of energy 
deposition in clusters. Conclusions will be drawn in the 
last section 



I. DESCRIPTION OF THE CODE. 

Let us start by describing the Classical N-Body Dy- 
namics code (labeled CNBD) used in this article. The 
basic ingredients of such a code are very simple. The dy- 
namical evolution of each particle of the system is driven 
by the classical Newtonian equations of motion. The 
two-body potential used in the present work is a third de- 
gree polynomial whose derivatives are null at the range 
n and at the distance of maximum depth r TO , n . The 
depth value is V m i n and the value at r = is finite and 
equal to Vq. This potential has the basic properties of 
the Lennard- Jones potential used in other works [TsUli^: 
a finite range attractive part and a repulsive short range 
part. To follow the dynamical evolution of the system 
an adaptative stepsize fourth-order Runge-Kutta algo- 
rithm is used [20j | . The main difference with other works 
is that the time step At can vary: if the potential varies 
strongly, At is small and when the potential varies gently, 
At becomes larger. This allows a very high accuracy with 
shorter CPU time than for fixed time step algorithms, ft 
requires an additional simulation parameter e which is 
adjusted to ensure the verification of conservation laws 
(energy, momentum, angular momentum) with a reason- 
able simulation time. The energy difference between the 
beginning and the ending simulation time is lower than 
0.001%. This simulation has five free parameters: four 
linked to the physics (the interaction) and one linked to 
the numerical algorithm. 

Since one wants to study the simplest case, neither long 
range repulsive interaction nor quantum corrections like 
a Pauli potential have been introduced [2lJ. Addition- 
ally, no statistical decay code is applied to the excited 
fragments formed during the collision. The final prod- 
ucts have to be regarded as "primary" products which 
would decay afterwards. Although most of the ingredi- 
ents necessary for a correct description of atomic nuclei 
are missing from this simulation, it should be noted that 
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the reaction mechanisms observed in nucleus-nucleus col- 
lisions are seen, and the properties of the "ground states" 
of such clusters are qualitatively close to those of nuclei 

In order to avoid any confusion with nuclear physics, 
the units used here are arbitrary and called Simulation 
Units (S.U.). The distance will then be in Distance Sim- 
ulation Units (D.S.U.), the energies in Energy Simula- 
tion Units (E.S.U.), the velocities in Velocity Simula- 
tion Units (D.S.U./T.S.U.) and the reaction time in Time 
Simulation Units (T.S.U.). 

For this study, 16,400 events have been generated. For 
a fixed projectile size N pro j, a fixed target size N targ and 
a fixed available energy in the center of mass E CAn jN, 
the impact parameter b is randomly chosen assuming 
a flat distribution between and b max — r(N pro j) + 
r(Ntarget) + where r(N pro j) is the mean square radius 
of the projectile, r(N tar get) is the mean square radius 
of the target and r\ the range of the potential. In the 
analyses, each event is weighted by its impact parameter 
value (weight oc b). The systems and the energies stud- 
ied here are summarized in tabled As it can be seen in 
table [I] the available energies in the center of mass frame 
Ec.m./N have been chosen in such a way that the excita- 
tion energy of the fused cluster E* /N is close to the bind- 
ing energy of the fused cluster (E c , m jN = 90 E.S.U.), 
close to the energy of the least bound particle in the fused 
cluster (E c . m ./N = 60 E.S.U.), far below (E c . m ./N = 
30E.S.U.) or far above (E c . m ./N = 120E.S.U.) the bind- 



ing energy of the fused cluster. The stars for the 100 + 
100 system mean that these energies were estimated from 
the liquid-drop parametrisation of the binding energies 
of the clusters [I3 ■ In that case, the energy of the least 
bound particle is taken equal to —65 E.S.U.. For each 
system, the excitation energy Ej used of the fused cluster 
is determined by: 



E*fused = Ec.m.+Ebind{N pro j)+Ei ) ind{Nta,rg)—Ebind(N fused) 

(1) 

where E^md {N) is the binding energy of the cluster of 
size N. 

In this unit system, the parameters of the reaction are 
the following: V = 540 E.S.U., V min = -20 E.S.U., 
"Train = 10 D.S.U. and 7"i = 15 D.S.U.. Typical mean 
square radii are around 10 D.S.U. for N « 20 and 
around 15 D.S.U. for N w 70 (for more information, 
see The time for a particle to go through a cluster 

ranges from 7 to 10 T.S.U. at the lowest available energies 
(Ec.m./N < 30 E.S.U/.) and ranges from 3 to 5 T.S.U. 
at the highest available energies (E c , m jN > 90E.S.U /.). 
As seen in fusion-like mechanisms and particle trans- 
fer mechanisms are observed at low E c . m jN values 
{E c . m ./N < 60 E.S.U.), while multi-fragment production 
and neck formation and break-up is dominant at high 
Ec. m ./N values (E c . m ./N > 60 E.S.U.). As it will be 
seen in section ILTll typical thermalisation times for such 
systems ranges from 15 to 20 T.S.U.. 
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FIG. 1: Excitation energy of clusters versus their parallel ve- 
locity for Nproj~13 on Ntarg=13 collisions at different avail- 
able energies in the center of mass. On each panel, the full 
line corresponds to the expected evolution for a pure binary 
process, the dashed horizontal line corresponds to the energy 
of the least bound particle for N — 26 and the circle corre- 
sponds to the expected values for the fused system. On each 
plot, the darkest gray regions correspond to the highest differ- 
ential cross section values. 



II. ENERGY DEPOSITION IN FINAL 
CLUSTERS. 

In this section, we will be interested in the energy de- 
position in clusters at the simulation time t = 2Q0T.S.U.. 
At this stage, the clusters are well separated from each 
other in configuration space. As it will be seen later, 
these clusters are thermally equilibrated. They can be 
viewed as primary clusters. For this first simple analy- 
sis, the clusters have been identified by using the Min- 
imum Spanning Tree method (labeled MST) which as- 
sumes that two particles belong to the same cluster if 
they are in potential interaction, i.e. if their relative dis- 
tance is below the range r\ of the potential. 

A. Variations of the energy deposition with the 
velocity of the cluster 



FIG. 2: Same as^\for N pro j=34 on Ntarg=34 collisions On 
each panel the dashed horizontal line corresponds to the energy 
of the least bound particle for N = 68. 
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On figures ^121 El EI and [S] are plotted the evolution of 
the excitation energy E* /N of the cluster with its par- 
allel velocity V/ / for the whole impact parameter range 
and for different available energies in the center of mass. 
The excitation energy of each cluster is simply the dif- 



FIG. 3: Same as^\for N pro j=18 on Nt arg =50 collisions On 
each panel the dashed horizontal line corresponds to the energy 
of the least bound particle for N = 68. 
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FIG. 4: Same as\^for N p 



-50 on N ta rg=50 collisions On 



each panel the dashed horizontal line corresponds to the energy 
of the least bound particle for N = 100. 



ference between the total energy (potential plus kinetic) 
and the "ground state" energy of the cluster: 



E * = j2 E? n + m *w*o (2) 

where E\™ is the kinetic energy of the particle i in the 
cluster's centre of mass, V(\rt — the potential energy 
between the particles i and j, Ebi nc [(N) the binding en- 
ergy of the cluster and N its number of particles. This 
excitation energy is determined at the end of the calcu- 
lation corresponding to t = 200T.S.U.. As it will be seen 
in the next section, this energy is very close to the one 
obtained at the separation time of the clusters (the small- 
est time at which clusters can be identified), since in this 
time range the emission of monomers or small clusters by 
the primary clusters (evaporation) is very weak and the 
clusters have no time to cool down significantly 01- O n 
each panel of the figures, the full line corresponds to the 
expected correlation between E*/N and Vn for a pure 
binary scenario, i.e. the excitation energy is only due to 
the velocity damping of each partner, without any par- 
ticle exchange between them. This excitation energy is 
then determined by: 
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FIG. 5: Same as^\for N pro j=100 on Ntar 3 =100 collisions 
On each panel the dashed horizontal line corresponds to the 
estimated energy of the least bound particle for N = 200. 



where m p and nit are the projectile and the target mass 
respectively, v re i_ max the initial relative velocity between 
the projectile and the target and v re i is the relative veloc- 
ity between the projectile and the target after interaction. 
The mass of each cluster is proportional to the number of 
particles m{N) = m par ticie x N where m par ucie = 20S.U.. 
The initial relative velocity is given by: 



m t 



m p m t 



(4) 



The small circle displayed on each panel is centered 
around the expected values of velocity and excitation en- 
ergy for the fused system. The horizontal dashed line 
corresponds to the energy of the least bound particle 
EheastBound for the fused system N = N proj + N target . 
This energy is determined for the stable clusters ( "ground 
state") and is defined as follows: 



E 



LeastBound 



max 



(5) 



where is the relative distance between the particles i 
and j and V(rjj) is the value of the two-body potential. 
As it has been shown in ^t|j this energy varies dramat- 
ically with N. In their "ground states" the clusters are 
small crystals and the EL east Bound value is mainly due 
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to geometrical effects (number of neighbors of a particle 
at the surface). Such variations are well known in cluster 
physics p^ . 

At low E c . m ./N values {E c . m ./N = 25 E.S.U. for the 
13+13 system and E c _ m _/N = 30 E.S.U. for the others), 
the points are slightly below the full line. This means 
that the excitation energy is strongly linked to the ve- 
locity damping. In that case, the collisions leads to the 
formation of excited projectile-like and target-like clus- 
ters. The small shift is due to mass transfers between 
the projectile and the target, and to promptly emitted 
clusters. The area corresponding to the fused system is 
well populated showing that a complete fusion process 
occurs. At intermediate energies {E c , m jN = 45 E.S.U. 
for the 13+13 system and E c . m ./N = 60 E.S.U. for 
the others), the distribution of points is roughly com- 
patible with the pure binary process hypothesis (forma- 
tion of excited projectile-like and target-like clusters) ex- 
cept for cluster velocities which lead to excitation ener- 
gies per particle higher than E LeastB ound{Nf use d) (where 
N f US ed = Np ro: j + N target ) in the pure binary process pic- 
ture. For these clusters, E* /N is always smaller than 
ELeastBound{N fused) ■ The complete fusion process area, 
located above E LeastB ound(Nf use d), is empty. In that 
case, as it will be seen in the next section, incomplete 
fusion process occurs. For the two highest energies, this 
trend is enhanced. Around the projectile and the target 
velocity the projectile-like and target-like clusters have 
an excitation energy compatible with the pure binary 
process hypothesis. Around the centre of mass veloc- 
ity, when this picture would give excitation energies per 
particle higher than E Least Baund(N fused), one finds clus- 
ters at small excitation energies. The energy of the least 
bound particle seems to be a limit to the excitation en- 
ergy which can be stored in these classical clusters. 



B. Variations of the energy deposition with the 
size of the cluster 



This limitation can be more clearly seen when the ex- 
citation energy E* /N is plotted as a function of N, as 
in figures El El M an d ED O n each panel, the full 
line corresponds to the energy of the least bound par- 
ticle ELeastBound in the cluster and the dashed line to 
the binding energy per particle E i n d/N . As in figures 
El El El an d El the small circle corresponds to the 
expected values for the fused system. At low energy, 
the area corresponding to complete fusion is filled and 
all the available energy can be stored as excitation en- 
ergy. But for higher energies, one can clearly see that for 
each fragment size, E* /N never overcomes Er, e astBound- 
At intermediate energy, clusters with sizes higher than 
the projectile size and the target size can be seen. This 
area corresponds to an incomplete fusion process. For 
the two highest energies, the plots are almost identical: 
there is no more fusion (the small circle area is empty) 
and the clusters are smaller than the target and than the 
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FIG. 6: Excitation energy of clusters versus N for N pro j —13 
on Ntarg=13 collisions at different available energies in the 
center of mass. On each panel, the full line corresponds to the 
energy of the least bound particle, the dashed and dotted line 
corresponds to the binding energy per particle and the circle 
corresponds to the expected values for the fused system. On 
each plot, the darkest gray regions correspond to the highest 
differential cross section values. 
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FIG. 7: Same asffifor N pro j =34 on Ntar S =34 collisions. 
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projectile. One can notice that E* /N never reaches the 
binding energy Eund/N except for small clusters where 
Ebind/N and E LeastBound are equal. The energy of the 
least bound particle in the cluster E^ eas tBound is an up- 
per limit for the energy per particle which can be stored 
in the cluster, whatever the system size, the available 
energy in the center of mass and the asymmetry of the 
entrance channel. 



C. Variations of the energy deposition with the 
impact parameter and the available energy in the 
center of mass. 

Figure ^] shows the variations of the average excita- 
tion energy per particle < E * /N > with the reduced 
impact parameter b re d — b/b max for different available 
energies in the center of mass frame E c m /N. The av- 
erage is calculated for clusters with TV greater or equal 
to 3. Each panel corresponds to a system. The squares 
correspond to the lowest E c . m ./N values, the triangles 
to Ec.m./N values close to ELeastBound, the diamond to 
Ec.m./N values close to Ebind/N and the circles to the 
highest Ec.m./N values. As in the previous figures, the 
average excitation energy is limited around AQE.S.U. ex- 
cept for the 100 + 100 system for which the maximum 
energy reached is 50 E.S.U. at E c . m ./N = 60 E.S.U.. 

At the lowest energy, for symmetric systems, < E * 
/N > increases when b re d decreases down to b re d = 0-4, 
and then is constant at < E * /N >« 40E.S.U. below. 
This saturation is due to the occurrence of fusion: the 



maximum excitation energy is reached by the fused sys- 
tem. For the asymmetric system 18 + 50, saturation 
occurs at smaller impact parameters. For intermediate 
energies, the picture is roughly the same. But although 
the available energy is higher, the maximum < E * /N > 
value is almost the same as for the lowest available en- 
ergy. In that case, an incomplete fusion process occurs 
and the excitation energy of the heaviest cluster is lim- 
ited by the energy of its least bound particle, which is 
well below the expected excitation energy of the fused 
system. For the two highest energies, < E * /N > in- 
creases when b re d decreases down to b re d — 0.4 and then 
< E * /N > decreases when b re d decreases. This effect 
results from the lower cluster sizes for the most central 
collisions, which at these energies produce several clus- 
ters of intermediate sizes. Indeed, the smaller is the clus- 
ter size, the lower is the upper limit in energy storage. 
Since clusters have a smaller size for central collisions, 
the < E * /N > decreases consequently. This evolution 
does not seem to depend strongly on the total system 
size or on the entrance channel asymmetry. 

The same evolution is seen on figure 1121 which shows 
the evolution of the ratio < E*/N > / (E* /N) ma x where 
(E * /N)max corresponds to the maximum excitation en- 
ergy expected for the fused systems. For low energies, 
this ratio continuously increases when b re d decreases and 
is close to 1 below b re d ~ 0.4: this corresponds to the oc- 
currence of the fusion process. This value of 1 is reached 
for b r ed below 0.2 for the 18 + 50 system. This difference 
can be understood quite easily: for a fixed Ec.m./N value, 
the relative velocity between the projectile and the target 
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FIG. 10: Same as\Qfor N pro j =100 on N ta rg=100 collisions. 
The dashed lines on each panel correspond to the extrapolated 
values of the binding energy (upper line) and the extrapolated 
energy of the least bound particle (lower line) for N > 100. 



is higher for an asymmetric system than for a symmet- 
ric one (see equation . This leads to an interaction 
time smaller for the asymmetric system than for a sym- 
metric one and then to less efficient energy deposition 
in the clusters. For higher energies, this value of 1 is 
never reached. The maximum value even decreases when 
Ec.m./N increases: the relative amount of energy stored 
in clusters is lower and lower when the available energy 
increases. 



D. Discussion 

These studies show that in classical systems, the en- 
ergy storage in clusters is limited. Such a saturation was 
observed for two-dimensional classical systems [2i| • This 
limitation is only linked to the intrinsic properties of the 
cluster: the energy deposition cannot be higher than the 
energy of the least bound particle in the cluster. 

This limitation of excitation energy can be understood 
quite easily. The least bound particle remains bound to 
the cluster only if its total energy is negative, i.e. its 
kinetic energy due to the excitation is below its potential 
one. If one assumes that the excitation energy is roughly 
equally shared over all particles in the cluster, when the 
kinetic energy balances the potential energy of the least 
bound particle, this particle is no longer bound to the 
cluster and quickly escapes. To be observed for a long 



time (greater than the thermalisation time), the excited 
cluster must have an excitation energy per particle below 
the energy of the least bound particle. 

The mechanism of energy deposition in classical N- 
body clusters seems to be the following one: the excita- 
tion is mainly driven by the velocity damping of the two 
partners and to a lesser extent by exchanges of particles 
between them. Once the energy of the least bound parti- 
cle is reached, unbound particles and/or clusters escape 
quickly and keep an excitation energy per particle below 
the energy of the least bound one. As a consequence, 
the highest energy deposition per particle can only be 
obtained at E c , m ,/N energies close to E LeastB ound- For 
higher available energies, the system fragments quickly, 
leaving rather "cold" clusters around the center of mass 
velocity. The energy in excess is stored in the kinetic en- 
ergies of the clusters. A similar mechanism could be an 
explanation of the apparent saturation of primary frag- 
ments excitation energies observed for central collisions 
of the Xe + Sn system when the incident energy increases 
from 32 to 50 A.MeV [3 El. This could also be an ex- 
planation for the limitation of temperature values around 
5 MeV by using isotope ratios methods or population ra- 
tios methods [jj, El ■ Such a limitation has been observed 
in the fragmentation of Uranium projectiles at relativis- 
tic energies This limitation was also suggested by the 
observation of the saturation of the evaporated neutron 
multiplicity when the incident energy increases 0]. 

Providing this conclusion can be applied in nuclear 
physics, i.e. if quantum effects and Coulomb interaction 
do not modify strongly the above picture, this means that 
the maximum energy per nucleon which can be stored 
in thermalised hot nuclear fragments would never exceed 
the energy of the least bound nucleon. This energy would 
correspond to the energy of the last populated level. If 
this assertion is true, that would be in contradiction 
with several experimental works in which high energy 
depositions (up to 20 A.MeV) were measured 0, EJ. 
This discrepancy may result from the calorimetry anal- 
yses used in these articles, whereas E * /N is directly 
extracted in this work. 

One has to be careful because the observations made 
here are for a fixed simulation time. At the early stages 
of the collision higher energy deposition could be reached 
for short times. The next section will be devoted to 
the evolution of this energy deposition with the reaction 
time. 



III. EVOLUTION OF THE ENERGY 
DEPOSITION WITH TIME. 

The aim of the present section is to check to which 
extent the limitation of the energy deposition observed 
for final (long simulation times) thermally equilibrated 
clusters is true during the collision. Since the limitation 
is linked to the energy of the least bound particle in the 
cluster, we will follow the ratio of the excitation energy 
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FIG. 11: Average excitation energy stored in the clusters as a function of the reduced impact parameter. The different lines 
and symbols correspond to different values of the available energy per particle in the center of mass frame. The E c . m ./N values 
indicated between parentheses correspond to the 13 + 13 system. See text for details. 



per particle to the energy of the least bound particle (E* 
/N)/E Least Bound(N) with collision time. 

Such an evolution is shown in figure IbH for MST clus- 
ters for central collisions (b re d < 0.1) for all systems at 
E c . m ./N = 90E.S.U. (65E.S.U. for the 13 + 13 system). 
It is seen that for MST clusters (open circles) very high 
energies can be found, up to the maximum excitation 
energy. This means that very excited systems could be 
formed at the early stages of the collision. But one has 
to check first if this transfered energy is equally shared 
between all degrees of freedom and hence if these inter- 
mediate clusters are thermalised. 

One way to check this thermalisation is to follow the 
ratio between the dispersion of the distribution of the ki- 
netic energy of particles o~(Eki n ) in the cluster and its 
average value < Ekm >• For a thermalised system, 
the value of this ratio is well defined: it is equal to 

\/2/3 in the canonical ensemble, and equal to \ 



in the micro-canonical ensemble [24], where N is the 
number of particles in the cluster. The time evolution 

of Rthermalisation — (^(^fein)/ ^ ^kin ^) / \J 3N— 1 ^ 

displayed on figure El f° r the 34+34 system (first row) 
and the 50+50 system (second row) at the two highest 
Ec.m./N values (90 E.S.U. left column and 120 E.S.U. 
right column) . A value close to 1 indicates that the clus- 
ter is thermalised. Let us first follow the evolution of 
the ratio for MST clusters (open circles). The ther- 
malisation of clusters is reached at t « 50 T.S.U. for 
both systems at E c . m ./N=90 E.S.U.. The excitation en- 
ergy per particle of the MST clusters at this time (see 
figure El open circles) is close to or below the energy 
ELeastBound of the least bound particle in the cluster. 
For later times, the ratio slightly varies around 1, show- 
ing that clusters stay thermalised even if evaporation of 
light clusters occurs. The same conclusion can be drawn 
at E c .m ,/iV=120 E.S.U. for which the thermalisation is 
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FIG. 12: Average amount of available energy stored in the clusters as a function of the reduced impact parameter. The 
different lines and symbols correspond to different values of the available energy per particle in the center of mass frame. The 
Ec.m./N values indicated between parentheses correspond to the 13 + 13 system. See text for details. 



reached at t ss 35 T.S.U. when the excitation energy per 
particle is very close to Er, eas tBound- From this simple 
study, one can conclude that the highly excited inter- 
mediate cluster is not thermalised, and that thermalised 
clusters have an excitation energy per particle close to or 

below ELeastBound- 

It has been found in earlier works that such systems are 
indeed composed of several small clusters at early times 
[Til l25l |2(| . The size distributions of these clusters have 
been found to be almost identical to those identified af- 
ter their separation [T9j . The problem is how to identify 
these clusters. In that respect, the MST algorithm is no 
longer suited since particles close to each other at these 
early stages are not necessarily bound at later times. Dif- 
ferent methods have been developed to identify clusters. 
One is to assume that two particles are belong to the 
same cluster if they are bound, i.e. if their relative energy 
is negative. These clusters have been labeled Coniglio- 



Klein clusters (CK clusters) or MSTE clusters in other 
works. Another way to define clusters is to find at each 
time step the most bound partition of particles in clus- 
ters. This algorithm is described in detail in j2?J and 
is called Early Cluster Recognition Algorithm (ECRA). 
The latter will be used in the present work. 

On figure U21 the evolution of (E * /N)/ E LeastBound 
with the collision time is shown for ECRA (filled circles) 
clusters. Above t w 50T.S.U. both algorithms (MST and 
ECRA) give the same results: clusters are well separated 
in the configuration space. At this time, clusters do not 
interact anymore with each other: this corresponds to 
the "freeze-out" configuration assumed in statistical mul- 
tifragmentation models [2^, EJ . One has to notice that 
this time coincides with the thermalisation time found 
in the previous paragraph. The thermalisation of ECRA 
clusters (see figure filled circles) is also reached at 
the "freeze-out" time. After this time, clusters have an 
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time (T.S.U.) 



FIG. 13: Variation of the ratio (E* /N)/ELeastBomd with time for various systems and for central collisions (b re d < 0.1) when 
the available energy in the centre of mass is close to the binding energy of the fused system. The open circles correspond to the 
MST clusters and the full circles to the ECRA clusters. 



excitation energy which is below EL east Bound (ratio be- 
low 1). But for earlier times, the configuration is more 
complex. While MST clusters reach very high excitation 
energies per particle, the ratio for ECRA clusters is al- 
ways around or below 1, except for the 100+100 system. 
Apart for the latter system, the limitation observed in 
the previous section seems to be still true all along the 
collision time. But what leads to the difference observed 
for the 100+100 system? 

Indeed, for such a heavy system, pre-formed clusters 
remain long enough together to interact strongly. If one 
of the particles escapes from one of these clusters, the 
probability that it is captured by another cluster is high. 
In that respect the clusters are acting as a confining wall 
and inhibiting quick escapes of the most energetic parti- 
cles. These interactions between pre-formed clusters pre- 
vent their thermalisation before the "freeze-out" time. It 
can be seen in figure 1131 that the maximum value of the 



ratio is higher for heavier systems. This effect should 
be reduced when the time during which clusters inter- 
act is reduced, i.e. when the available energy increases. 
This is observed in figure 1151 which shows the evolu- 
tion of the ratio with time for E c . m jN — 120 E.S.U. 
(E c . m ./N = 85 E.S.U. for the 13 + 13 system). At this 
energy, the reaction times are smaller and the ratio is 
always close to or below 1. 

It could be argued that the observed evolution of 
ECRA clusters is only due to the algorithm of cluster 
recognition used. Since one tries to find the most bound 
partition, this could minimize the energies of clusters. 
It is indeed observed that ECRA clusters have an exci- 
tation energy below the MST clusters. One can make 
two remarks. The first one is that the EL east Bound has 
not been used to determine the ECRA clusters. There 
is hence no reason to find -ELeastBound as an upper limit 
for the excitation energy for ECRA clusters. The sec- 
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FIG. 14: Variation of the ratio (a(Ekin)/ < Ekin >)/y 3JV-2 time for the 34 +34 system (first raw) and the 50+50 system 

(second raw), for central collisions (bred < 0.1) at the two highest E c . m ./N values (90 E.S.U. left column and 120 E.S.U. right 
column). The open circles correspond to the MST clusters and the full circles to the ECRA clusters. 



ond remark is that ECRA clusters with energies higher 
than E Least Bound can be found in some specific cases (see 
figure IT31 for the 100 + 100 system) when particles can- 
not easily escape from clusters. The evolutions seen for 
ECRA clusters are then more likely due to the physics of 
the collision rather than to the definition of these clus- 
ters. 

As in section m the energy which can be stored in a 
thermalised cluster seems to be limited almost through- 
out the collision time. The detailed analysis shows that 
the exact value of this limitation may depend on reac- 
tion time and on the system size. Before the "freeze- 
out" time, this limit can be higher than ELeastBound if 
clusters are close enough together for a long time. This 



has been observed for the heaviest systems and for an 
Ec.m./N value close to the binding energy of the fused 
system. After the "freeze-out" time, all clusters have an 
excitation energy per particle below ELeastBound- 



If clusters are excited in a confining "box" (wall, 
neighboring clusters) which prevents quick emission of 
monomers and small clusters, very high excitation en- 
ergies can be deposited in these clusters: they are ar- 
tificially bound by the confining wall or the neighbor- 
ing clusters. If clusters are "free" (no confining wall, no 
neighboring clusters), the energy deposition per particle 
can not exceed E Le astBound- 
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FIG. 15: Variation of the ratio (E* /N)/ELeastBomd with time for various systems and for central collisions (b re d < 0.1) when 
the available energy in the centre of mass is close to the most bound particle in the fused system. The open circles correspond 
to the MST clusters and the full circles to the ECRA clusters. 



IV. CONCLUSIONS 



The energy deposition in clusters in cluster-cluster col- 
lisions has been studied in the framework of Classical 
N Body Dynamics. The mechanism of energy deposi- 
tion seems to be the following: the excitation energy of 
clusters is mainly due to the relative velocity damping 
between the projectile and the target, and to a lesser 
extent to particle exchanges between them. The ex- 
citation per particle of thermalised clusters is limited 
by the energy of the least bound particle in the clus- 
ter. In the low energy regime {E c . m JN < E Least Bound), 
the excitation energy increases when the impact param- 
eter decreases and reaches its maximum value for the 
most central collisions (b re d < 0.5). For the intermedi- 
ate energy regime (E c , m ,/N w E LeastBo und) the picture 
is almost the same except that for central collisions the 
incomplete fusion system has an excitation energy per 



particle close to E LeastBo und- For high energy regimes 
{Ec.m./N > Esind/N), the maximum energy deposition 
is found for intermediate impact parameters (b re d « 0.5) 
whereas for central collisions clusters are less excited and 
E * /N is always lower than E Least Bound- 

This limitation of energy deposition is almost verified 
throughout the collision. High energy deposits (E* /N w 
1.4 EL eas tBound) have been found for ECRA clusters and 
for heavy systems before the "freeze-out" time, which 
coincides with the thermalisation time of the clusters. 
This effect is mainly due to the times scales of the reac- 
tion. Above the "freeze-out" (thermalisation) time the 
excitation energy per particle of free clusters is below 

E L east Bound- 

This limitation of energy deposition in thermalised 
clusters could be an explanation for rather low excitation 
energies found in primary fragments in nucleus-nucleus 
collisions. If this limitation is also true for nuclei, this 
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would be in contradiction with many experimental works tions. 
in which excitation energies well above the binding en- 
ergy have been determined. This could shed new light on 
caloric curve analyses and on phase transitions studies. 
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